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ABSTRACT 

We present new ~ 1" resolution data of the dense molecular gas in the central 50 — 
100 pc of four nearby Seyfert galaxies. PdBI observations of HCN and, in 2 of the 4 
sources, simultaneously HCO + allow us to carefully constrain the dynamical state of 
the dense gas surrounding the AGN. Analysis of the kinematics shows large line widths 
of 100-200 km/s FWHM that can only partially arise from beam smearing of the 
velocity gradient. The observed morphological and kinematic parameters (dimensions, 
major axis position angle, red and blue channel separation, and integrated line width) 
are well reproduced by a thick disk, where the emitting dense gas has a large intrinsic 
dispersion (20-40 km/s), implying that it exists at significant scale heights (25-30% 
of the disk radius). To put the observed kinematics in the context of the star burst 
and AGN evolution, we estimate the Toomre Q parameter. We find this is always 
greater than the critical value, i.e. Q is above the limit such that the gas is stable 
against rapid star formation. This is supported by the lack of direct evidence, in these 
4 Seyfert galaxies, for on-going star formation close around the AGN. Instead, any 
current star formation tends to be located in a circumnuclear ring. We conclude that 
the physical conditions are indeed not suited to star formation within the central 
~ 100 pc. 

Key words: galaxies: kinematics and dynamics - galaxies: nuclei - galaxies: Seyfert 
- galaxies: ISM - galaxies: individual (NGC 2273, NGC 3227, NGC 4051, NGC 6951) 



1 INTRODUCTION 

The key component of the unification model for active 
galactic nuclei (AGN) is a geomet r ically thick torus (e.g. 
lAntonuccil 1 19931 . lUrrv fc Padovanil Il995l ) where gas and 
dust obscure the primary optical-ultraviolet AGN emission, 
depending on the line-of-sight viewing angle. The signature 
of such a torus is a pronounced i nfrared (IR) pea k in the 
spectral energy distribution (e.g. lElvis et al.lll994l ). which 
is interpreted as thermal emission from hot dust heated 
by the primary radiation (|Barvainidli987l ). In this scenario 
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the obscuring material surrounding the AGN accretion 
disk has a hydrostatic toroidal geometry and is thought 
to begin at the dust sublimation radius a nd extend for 
tens, or possibly even hundred s , of parsecs (iPier fc Krolikl 
1 19921 : [ Granato fc Danesd Il994l: INenkova. Ivezic. fc Elitzu 



2002: iNenkova et alj 120081 : ISchartmann et al.l 
iHonig et al.l 120061 ). A smooth continuous 



20051 . 1200S 
distribution 



of the medium in the torus leads to stability prob- 
lems. Thus it has been sugge sted that the thick torus 
shoul d either be clumpy (Nenkova. Ivezic. fc Elitzurl 



2002 



lElitzur. Nenkova, fc Ivezid 



2004; 



Fritz. Franceschini. fc Hatziminaogloul 120061 ). or be sup- 



ported by an additional force combined with the thermal 
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pressure. Nuclear star formation (SF) can play an important 
role in sustaining i ts vertical extent, e.g. through stellar 
radiation pressu re l|Thompson. Quataert^ fc Murray! 120051; 
Ballan tvne 2008]), super novae explosions JWada fc Norman! 



20021 ). or stellar winds |Navakshin fc Cuadral 12007! ). Thus 



the connection between AGN activity and nuclear SF is an 
important issue related to the torus s tructure. Indeed nu- 
clear SF happens in all types of AGN (jCid Fernandes et al.l 
12004 ) . and there is evidence that enhanced SF (reaching 
starb urst intensities) is related to the black hole accretion 
rate l|Sani et al.ll2010l ). 

Thanks to the capabilities of the adaptive optics as- 
sisted integral field spectroscopy, it has recently become pos- 
sible to study the in ner 100 pc o f galax ies with a spatial 
resolution of ~10 pc. lDavies et al.l (|2006h used SINFONI to 
analyse the molecular gas and stars in NGC 3227, finding 
the central gas distribution was geometrically and optically 
thick, analogous to the standard molecular torus. Further- 
more they found a recent, but no longer active, starburst 
within the tens of parsecs around the AGN that occurred 
~ 40 Myr ago. A similar result is found in a larger sam- 
ple of AGN galaxies, where there appears to be a delay of 
50-100 Myr between the onset of star for mation and the on- 
set of AGN activi ty ([Da vics ct al. 2007]). On similar scales 
of tens of parsecs, iHicks et al.l l|2009l . hereafter H09) traced 
the distribution and kinematics of the warm molecular gas, 
which they associated with the largest scale structures of the 
torus. Analysis of the H2 (1-0) S(l) emission in these AGN 
showed that the gas has a huge velocity dispersion (50-100 
km/s), implying it has a geometrically thick structure with 
a height of the order of the disk radius (H ~ R). Moreover, 
based on spatial distribution, column density and kinemat- 
ics, H09 concluded that the molecular gas is spatially mixed 
with the nuclear stellar population. 

In an attempt to further understand the circumnu- 
clear structure around AGN, we present here HCN(l-O) 
and HCO + (1-0) distributions and kinematics in four nearby 
Seyfert galaxies observed with the Plateau de Bure Interfer- 
ometer (PdBlfl With this array, it is possible to reach a 
spatial resolution of ^ 1" at 89 GHz. Although this is still 
one order of magnitude larger than the resolution obtain- 
able at near-IR wavelengths, the millimetre regime offers a 
number of clear advantages, most importantly the excellent 
spectral resolution and that the emission is free of extinc- 
tion. 

Our main purpose of this paper is to understand 
whether the high dispersion seen previously in the warm 
H2 is related to the tracer properties, or whether it is 
really representative of the bulk gas properties. This is an 
important distinction because the 2.12 /im H2 line traces 
just the warm (~ 1000 K) component of the gas, which 
accounts only for a tiny fraction of the total gas mass 



(10- 



10 -5 ). In contrast, the HCN(l-O) and HCO+(1-0) 



transitions probe the cold (~ 10 — 100K) and dense 
(n > 3 x 10 4 cm~ 3 ) mol ecular gas jGracia-Carpio et al.l 
l200rj ; iPapadopouiosl 120071 ) and are thus suitable for our 
purposes. This is the case mainly for two reasons: (i) 



HCN and HCO+ lines trace ~ 100 - 500 times denser 
gas than corresponding (rotational) CO transitions and 
(it) HCN (and HCO + ) are relatively strong close to 
the AGN, and weaker elsewhere. In fact the HCN (and 
HCO + ) to CO intensity ratio can be significant ly higher in 
AGN t han in starburst or quiescent regions (IKrips et al.l 
20071: [Gracia-Carpio. Garcfa-Burillo, fc Planesas! 120081 : 



Kohno et all I2OO8I : IDavies, Mark, fc Sternberg) |2012| ). We 
note that the X-ray irradiation of the gas by the AGN may 
affect the nuclear chemist ry (i.e. molecular a bunda nces). 
For example, mode ls in Boger fc Sternberg] |2005j) and 
IMeiierink fc SpaansI l|2005 ) show that the equilibrium 
abundances depend on the ratio of the local gas density and 
the incident X-ray ionization rate. These effects are beyond 
the scope of this work, but ar e investigated with specific 
refere nce to NGC 3227 by IDavies. Mark, fc Sternberg! 
|2012l ). Despite this, a modified chemistry should have little 
impact on the molecular gas kinematics. The use of two 
tracers, when available, reduces the risk that the observed 
kinematics are skewed by chemistry affecting the inferred 
gas mass distribution. As such, HCN and HCO + lines can 
be used to reliably measure the characteristic kinematics of 
the dense gas in the central region. 

In this paper we present HCN(l-O) observations for 
NGC 2273, NGC 3227, NGC 4051 and NGC 6951 together 
with HCO + (1-0) data for NGC 2273 and NGC 4051. Ober- 
vations and data reduction are discussed in Section 2, and 
the general properties of the molecular gas are presented in 
Section 3. In Section 4 we analyse the dense gas kinematics, 
and Section 5 is dedicated to the evolutionary interpreta- 
tion. Our conclusions are given in Section 6. 



2 OBSERVATIONS 

We obtained data using the PdBI (which has 6 antennaes of 
15 m diameter) in the extended A configuration (760 m max- 
imum baseline) for two sets of observations, the first relating 
to NGC 3227 and the second for NGC 2273 and NGC 4051. 
NGC 6951 was observed previously in A+B configuration, 
and the data published by K07. However, the K07 work 
deals mostly with the spatial trend of the HCN/HCO + ra- 
tio. Our analysis here concentrates instead on the central 
source kinematics. All the data presented here were pro- 
cessed and calibrated using the CLIC program in the IRAM 
GILDAS package^. We refer the reader to K07 for all the 
details concerning the observing conditions for NGC 6951. 
The observations of the other targets are described below. 

Observations of NGC 3227 were carried out during Jan- 
uary and February 2009 for program S098. In the 1 mm 
band, the H 13 CN(3-2) line at 259.01 GHz was observed in a 
single 7.6 hr track using the dual polarization receivers. At- 
mospheric conditions were moderate on this night and the 
system temperature (T 3y s) was about 160 K. In the 3 mm 
band the H 12 CN(l-0) line at 88.63 GHz and the H 13 CN(l-0) 
line at 86.34 GHz were observed simultaneously. This was 
achieved by observing a 1 GHz bandwidth segment across 
each line using a single polarization for each of the two 



1 The PdBI is managed by the IRAM institute, which is sup- 
ported by INSU/CNRS (France), MPG (Germany) and IGN 
(Spain). 



2 The GILDAS software is available at 
http: / /www.iram. fr/IRA MFR/GI LDAS] 
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Source Program Beam Size Beam PA Channel res. rms pc scale 

" X " deg km/s mjy/beam pc/" 



NGC 2273 


T095 


1.19 x 0.90 


38 


21 


0.76 


130 


NGC 3227 


S098 


1.19 x 0.68 


47 


50 


0.65 


85 


NGC 4051 


T095 


1.07 x 0.72 


70 


21 


0.63 


50 


NGC 6951 


PB67 


1.42 x 1.09 


73 


4 


1.6 


102 



Table X. (1) Source name. (2) Observing program. (3) Beam size in A configuration, or A+B configuration for NGC 6951. (4) Position 
angle of the beam, measured East of North. (5) Spectral resolution (after binning). (6) Root mean square uncertainty. (7) Spatial scale 
at the distance of the source. 



Source 


Flux Density 


Major x Minor axis 


Pos. angle 




mjy 


" X " 


deg 


NGC 2273 


2.34 ±0.09 


(1.86±0.05) x (1.13±0.03) 


87±3 


NGC 3227 


1.5 ±0.3 


(1.15±0.06) x (1.05±0.08) 


47±7 


NGC 4051 


0.9 ±0.2 


(1.13±0.07) x (0.90±0.09) 


82±10 



Table 2. Parameters of the continuum emission for the new observations. (1) Source name. (2) Flux density. (3) Observed major and 
minor axis diameters. (4) Position angle of the major axis, measured East of North. 



segments. These data were acquired over tracks totalling 
14.5 hrs during 3 nights with mediocre and variable atmo- 
spheric conditions, with up to 4-6 mm water vapour and 
sometimes strong wind. The T sy s was 80-100 K. The band- 
pass of the receivers was calibrated by observing the quasar 
3C 273 (first and second night) and 0851±202 (third night). 
Amplitude and phase calibrations were achieved by monitor- 
ing 0923±392, whose flux densities were determined relative 
to 3C 273, MWC349 or 0851±202 resepctively. We applied 
a spectral binning of 50 km/s to increase the signal to noise. 
The resulting angular beam sizes were 1.19" x 0.90" at posi- 
tion angle PA= 38° (East of North) for the H 12 CN(l-0) line; 
1.21" x 0.93" at 37° for the H 13 CN(l-0) line; 0.56" x 0.38" 
at PA= 23° for the H 13 CN(3-2). 

A detection of the H 1 2 CN(l-0) line in NGC 3227 
was pr eviously reported by ISchinnerer. Eckart. fc Tacconil 
(2000a) at a resolution of 2.4". Our more sensitive and 
higher resolution observations resolve the source and de- 
tect the 3 mm continuum. Neither the H 13 CN(l-0) nor the 
H 13 CN(3-2) lines were detected. 

We obtained 3 mm data for the Seyfert nuclei NGC 2273 
and NGC 4051 with program T095, executed during Decem- 
ber 2009. The observations were performed as for NGC 3227 
but with the central frequency shifted in order to simultane- 
ously observe both HCN(l-O) at 88.63 GHz and HCO + (1-0) 
at 89.19 GHz. NGC 2273 was observed with 2 tracks of 1.5 
+ 5.5 hr on source. For both tracks, the atmosphere was sta- 
ble, resulting in T sys = 65 K and T 3y3 = 80 K respectively. 
A single track of 7.7 hr was obtained for NGC 4051 with 
good sky conditions (wind > 15 m/s just at the end of the 
track) and T 3ys = 50 K. For NGC 2273, the bandpass was 
calibrated by observations of 0923±392 (short track) and 
3C454.3 (long track), while 0646±600 and 0716±614 were 
used as amplitude and phase calibrators respectively for the 
two tracks. During the long track, fluxes were calibrated 
using 3C454.3 and checked on 0716±714 and 0646+600; 
while for the short track fluxes where calibrated by mon- 
itoring 0646±600, 0716±714 and 1055±018 with respect 
to 0923±392. For NGC 4051, the bandpass calibrator was 



3C273. Phase and amplitude calibrations were achieved by 
monitoring J1224±435 and 1128±385, whose flux densities 
were determined relative to 0923±392. The data quality as- 
sessment allowed us to apply a 20 km/s spectral binning. 
The final angular beams were: 1.19" x 0.68" at PA= 47° 
and 1.07" x 0.72" at PA= 70° in NGC 2273 and NGC 4051 
respectively. The rms phase noises were between 30° and 43° 
at 3 mm, which introduces absolute position errors 0.1" 
in our observations. But we note that atmospherically or 
instrumentally induced relative position errors within each 
dataset are significantly smaller than this, so that separa- 
tion measurements depend instead almost entirely on the 
centroiding accuracy. 

To achieve a high spatial resolution, here we analyse 
uniformly weighted maps. This allows us to avoid contam- 
ination by circum nuclear structures, such as the ring emis- 
sion in NGC 3227 (|Davies et alj|2006r i and NGC 6951 (K07). 
Table 1 lists the observing parameters. 



3 RESULTS 

We first describe the overall properties of our sources, fol- 
lowed by a more detailed analysis source by source. Fig- 
ures[T][3] and[7] show the continuum, HCN, and (when avail- 
able) HCO+ maps for NGC 2273, NGC 4051 and NGC 3227 
respectively. For these sources, because the continuum is 
only weakly detected in single spectral channels, we obtain 
continuum maps by collapsing all the channels with no sig- 
nal from the molecular transitions (i.e. neither HCN nor 
HCO + ). The line maps are created by summing those chan- 
nels where the emission is detected and subtracting the con- 
tinuum. In Fig. [9] we show the HCN emission for NGC 6951 
from K07. 

The 3 mm continuum is detected in 3 of the 4 objects 
as a single compact (mostly) unresolved source. The sizes 
and flux densities of the continuum sources are given in Ta- 
ble Similarly, the line emission appears with a centrally 
concentrated morphology located at the same position, but 
(marginally) resolved. We measure the extent of each source 
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Molecular Lines 



Source 


Line 


Flux 

Jy km/s 


Major x Minor axis 
" x " 


PA 

deg 


r/b Sep 
// 


PA r/ b 
deg 


FWHM 

km/s 


NGC 2273 


HCN 
HCO+ 


1.76 ±0.07 
2.13 ±0.08 


(1.87±0.05) x (1.17±0.03) 
(1.94±0.04) x (1.34±0.04) 


18±6 
44±6 


1.10±0.02 
0.87±0.03 


23.2±0.4 
26.2±0.8 


181 ± 21 
175 ± 16 


NGC 3227 


HCN 


1.86 ±0.27 


(1.28±0.06) x (1.04±0.06) 


47±2 


0.38±0.01 


-37.1±0.4 


207 ± 34 


NGC 4051 


HCN 
HCO+ 


0.91 ±0.05 
2.07 ±0.05 


(1.92±0.09) x (1.02±0.05) 
(1.65±0.05) x (1.27±0.04) 


80±5 
72±2 


0.55±0.02 
0.60±0.04 


39.4±0.6 
19.6±0.7 


74 ± 10 
90 ± 10 


NGC 6951 


HCN 


1.02 ±0.02 


(1.44±0.03) x (1.10±0.03) 


-80±3 


0.60±0.01 


-50.2±0.5 


173 ± 29 



Table 3. Parameters of the molecular emission: (1) source name. (2) line transition: HCN(l-O) or HCO + (1-0). (3) Line flux within the 
central 3". (4) Observed major and minor axis diameters of the molecular emission. (5) Position angle of the major axis, measured East 
of North. (6) Separation of the red and blue velocity channel maps (see text for details of the channels summed to produce these). (7) 
Position angle (east of north) of the vector joining the centroids of the red and blue channel maps. (8) Width of the molecular line 
(FWHM) from the integrated spectrum. 



NGC 2273 . 
3mm continuum 



■ 



■A 




Figure 1. The 3 mm continuum (left, rms = 0.04 Jy beam km s -1 ), HCN (middle, rms = 0.08 Jy beam km s — 1 ) and HCO + (right, 
rms = 0.11 Jy beam - 'km s — 1 ) emission maps of NGC 2273. The contour levels are at 2, 4, and 6, times their respective noise levels. The 
beam size (FWHM) is shown in the lower left of each panel, and the white 'plus' signs denote the peak of the continuum. The continuum 
appears partially resolved, while the line emission is more obviously extended in both HCN and HCO + lines. 



both as the projected shape in the uv plane and in spatial 
coordinates. Because the data are characterized by a rela- 
tively low signal-to-noise ratio (S/N), the fitting in uv coor- 
dinates is performed using a symmetric Gaussian. This bi- 
dimensional fitting is performed using the MAPPING pro- 
gram in the GILDAS package, and gives a direct measure 
(i.e. without beam convolution) of the intrinsic source size. 
If the source distribution is highly non-circular this may not 
be realistic and so this size estimate is used only as a consis- 
tency check. Our primary size measurement is based on a fit 
to the data in spatial coordinates (i.e. including the cleaned 
beam). Thus, the FWHMs in the two spatial directions rep- 
resent the observed major and minor axes of the sources. We 
note, however, that the fluxes and their relative uncertain- 
ties are integrated within a 3" aperture around the image 
center. To properly compare the observed source dimensions 
(in spatial coordinates) with intrinsic measurements (in the 
uv plane), it is necessary to subtract in quadrature the beam 
shape from the former. We derive the uncertainties for the 
observed sizes using Monte Carlo techniques by creating 10 4 
realizations of the image, each time adjusting the data val- 
ues with random numbers according to the measured noise 
and re-fitting the Gaussian and background. For the two ob- 



jects in which both lines are observed, the similarity of the 
spatial parameters of the HCN and HCO + emitting regions 
given in Table [3J verifies that the two molecules are likely to 
be tracing the same gas components. 

The spectra of the sources are shown in Figures [2] [4l [8] 
and llOl These are integrated over all spatial elements within 
a 3" aperture in which the signal exceeds a 0.25 mjy beam - 
threshold. This threshold is significantly lower than the 
rms given in Table 1 and is simply used to increase the 
S/N without affecting the resulting integrated line widths 
(which are consistent within the errors, both with and with- 
out the threshold). We derive the uncertainties for the line 
FWHMs using Monte Carlo realizations of each spectrum, 
as described above for observed source sizes. The most re- 
markable characteristic is the large line width: 207 ±34 and 
173±29 km/s for the HCN line in NGC 3227 and NGC 6951, 
respectively; and mean FWHMs of 177±14 and 82±10 km/s 
for NGC 2273 and NGC 4051 respectively. Such widths are 
comparable with those observed for the warm gas traced by 
the H2 (1-0) S(l) line, reported by H09. As discussed in the 
next section, the implied high velocity dispersion is a crucial 
parameter to probe the dense gas distribution. 

The spatial sizes, velocity widths, and fluxes of the lines 
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HCN(1 -0) 




-1000 -500 

Velocity offset (km/s) 



HC0+(1 -0) 




-500 500 1000 1500 

Velocity offset (km/s) 

Figure 2. Integrated spectra of NGC 2273 for the HCN (top), 
and HCO + (bottom) lines. The blue curve is a fit (Gaussian plus 
constant term) to the 3 mm continuum and the line emission, 
which can be seen to have a large velocity width (mean FWHM 
of the two lines is 177 ± 14 km/s). The dashed line represents the 
zero flux level. 



are listed in Table [3] We emphasize that the measurements 
in this section are observed quantities. To perform an appro- 
priate dynamical analysis of the gaseous structure we need 
intrinsic properties, i.e. we have to take into account the 
impact of spatial and spectral beam smearing (see Section 4 
for details). 

3.1 NGC 2273 

The 3 mm continuum appears to be partially resolved, hav- 
ing a FWHM of 1.86" x 1.13" at a PA= 87° (from a Gaussian 
fit to the data in the left panel of Fig. [TJ. The 3 mm flux 
density of 2.46 ± 0.16 mjy from a Gaussian fit to the con- 
tinuum map is consistent with the mean spectral value of 
2.34 ± 0.09 mjy from the integrated spectrum (Fig. [2]). 

The HCN line (Fig. [1] middle panel) has an integrated 
flux of 1.76 ± 0.07 Jy km/s in a 3" aperture. Its FWHM is 
1.87" x 1.17" at a position angle of 18°. Allowing for the 
beam size, this is consistent with the projected intrinsic size 
of 1.3"±0.3" from a symmetric Gaussian fit to the uv points. 
Thus the line emitting region is compact but nevertheless 
spatially resolved. 

This conclusion is confirmed by comparing the HCN 
properties with those of the HCO + line (Fig.[T] right panel), 
which has an integrated flux of 2.13 ± 0.08 Jy km/s in a 
3" aperture. The latter has a FWHM of 1.94" x 1.34" that, 
given the errors derived with Monte Carlo realizations listed 
in Table O is consistent with the corresponding measure- 




500 1000 

Velocity offset (km/s) 



Figure 4. Integrated spectrum of NGC 3227 showing the 3 mm 
continuum and the HCN line with a Gaussian fit (in blue) showing 
the large velocity dispersion of the line (FWHM=207 km/s). 



ment above for the HCN line with respect to the major 
axes, but implies a more extended source along the minor 
axis. A symmetric Gaussian fitted to the uv table gives an 
intrinsic FWHM of 1.7" ± 0.4", also in agreement with the 
size based on spatial coordinates. 

The position angle of the HCO + major axis appears to 
differ slightly from that of the HCN, and their centers are 
marginally offset with respect to the continuum. This could 
reflect different distributions of the tracers even though the 
offset observed is of the same order of the positional accu- 
racy. On the other hand, the PAs of the red/blue channels 
for the two lines (see Sec. [3} are in much better agreement. 
The HCN and HCO + lines are also similar spectrally. Indeed 
a single Gaussian fit to the spectra in Fig. Ogives FWHMs 
of 181 km/s and 175 km/s for HCN (top panel) and HCO+ 
(bottom panel) respectively. We have derived uncertainties 
of 21 km/s and 16 km/s with Monte Carlo techniques. 



3.2 NGC 3227 

A Gaussian fit to the 3 mm continuum (Fig [3] left panel) 
yields a size of 1.15" x 1.05" FWHM at a position angle 
of 47°. This is very similar to the beam size and indicates 
the continuum is spatially unresolved. The 3 mm flux den- 
sity from the integrated spectrum (Fig2} is 1.56 ±0.24 mJy, 
consistent with the 1.79T0.12 mjy obtained from the Gaus- 
sian fit to the continuum map. 

The HCN line has an integrated flux of 1.86 ± 
0.27 Jy km/s in a 3" aperture (Fig [3] right panel). Com- 
parison with th e previous measurement of 2.1 Jy km /s at 
2.4" resolution (|Schinnerer. Eckart. fc Tacconil [20~00at ) sug- 
gests that very little of the line emission has been resolved 
out at our higher resolution, and that most of the HCN in 
NGC 3227 does originate from the central compact source. 
The compact nature of the HCN emission is in stark con- 
trast to the CO(2-l) emission, as shown in Fig. [5] Indeed 
the CO emission is distributed around the circumnuclear 
ring (whi ch is also seen in the H-band stellar continuum, 
Fig 5 and lDavies et al]|2006l ). and very little originates from 
the central arcsecond. On the other hand, the HCN emission 
is dominated by the nucleus itself. 

A symmetric Gaussian fit to the data in the uv plane 
yields a projected intrinsic FWHM of 0.8"±0.3". This is con- 
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1 

1 1 

1" 


NGC 3227 ' 
3mm continuum 




HCN(l-O)' 




_ 






Is 
















■ 


1 


N 

E J 

1 




1 




i 



Figure 3. The 3 mm continuum (left, rms = 0.07 Jy beam _1 km s" 1 ) and HCN (right, rms = 0.11 Jy beam 'km s 1 ) emission maps 
of NGC 3227. Labels are as in Fig. [T] The contour levels in the left panel are at 3, 5, and 7 times the noise level; in the right panel at 2, 
4, and 6 times the noise level. The white 'plus' sign in the right panel indicates the location of the continuum peak. The continuum is 
unresolved while the line emission is clearly extended. 




RA offset (arcsec) 

Figure 5. Maps of the central few arcsecs of NGC 3227 (1" ~ 80 pc), adapted from iDavies et al] l|2006t l. In each panel the dashed 
ellipse traces the circumnuclcar ring, and the square indicates the field of view of the SINFONI data presented by these authors. Left: 
H-band image of the circumnuclear ring (obtained after su btracting a bulge and disk model fitted to larger scales). Center: the CO(2-l) 
emission (data from Schinncrcr, Eckart, & Tacconi 2000a) is distributed around the ring, with relatively little from the nucleus. Right: 
the HCN(l-O) peaks at the nucleus with relatively little originating in the ring. North is up and East is to the left (i.e., negative offsets). 



sistent with the FWHM of 1.5" x 1.04" (at PA 47°) measured 
from the reconstructed image, once the finite beam size is 
taken into account. It indicates that the line emission is re- 
solved. However, we note that the long axis of the nuclear 
emitting region coincides with the minor axis of the circum- 
nuclear ring, which as Fig. [S]shows is traced by the CO(2-l) 
emission. Along this axis, HCN emission from the ring may 
be blended with that from the nucleus, which would bias 
the size measurement of the nuclear source. We have there- 
fore separated these contributions by extracting a profile 
along this position angle (see Fig(6| and fitting 3 indepen- 
dent Gaussians at fixed positions - representing the nucleus 
and a cut through the ring on either side - to the overall 
HCN profile. One of these, which would be associated with 
the ring to the south-west, has a negligible contribution and 
so does not appear in the plot. However, the north-east side 
of the HCN profile is clearly broadened by a subsidiary peak. 
Since the asymmetry of the full profile can be matched by 



the addition of a component at the same location as the CO 
(which arises from the ring), we conclude that it is associ- 
ated with the ring. Accounting for this reduces the observed 
FWHM of the nuclear component slightly to 1.28". Quadra- 
ture correcting the observed 1.28" x 1.04" size of this com- 
ponent for the beam indicates that the intrinsic source size 
may be as small as 0.5" along both axes, a little less than 
implied by direct fitting of the uv data. 



As found in NGC 2273, the HCN line in NGC 3227 
has a remarkably large velocity width. A simple Gaussian 
fit to the integrated spectrum in Fig. 2] gives a FWHM of 
207 ± 34 km/s, where, as before, we have used Monte Carlo 
techniques to estimate the uncertainty. 
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Figure 7. The 3 mm continuum (left, rms = 0.02 Jy beam - x km s — 1 ), HCN (middle, rms = 0.04 Jy beam - x km s — 1 ) and HCO + (right, 
rms = 0.05 Jy beam - 1 km s — x ) emission maps of NGC 4051. Labels are as in Fig.fTJ The contour levels run between 2—6 times (left) and 
2-8 times (middle) the noise level in steps of 2 <r; and between 3-12 times (right) the noise level in steps of 3 a. The white 'plus' signs 
denote the location of the continuum peak. The continuum is basically unresolved while the line emission is clearly extended. 
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Figure 6. NGC 3227 profiles of the CO(2-l) (solid blue line) and 
the HCN(l-O) (solid red line) extracted in a 0.6" wide slit along 
a position angle of 45° (i.e. along the minor axis of the circum- 
nuclear ring). The CO peaks in the ring, while the HCN peaks 
on the nucleus. Superimposed is a fit to the HCN profile (dashed 
green line) which comprises 3 Gaussians (dotted black lines; one 
of the Gaussians is not seen because it has zero flux). The sec- 
ondary peak on the right, associated with the ring, matches the 
broad asymmetric overall profile. Taking just the central Gaus- 
sian yields an observable size for the nuclear emission of 1.28" at 
this PA. 

3.3 NGC 4051 

A Gaussian fit to the continuum map of NGC 4051 (Fig. [7] 
left panel), gives a FWHM of 1.13" x 0.90" at PA= 82°. 
Comparing this to the 1.07" x 0.72" beam at 80° suggests 
the 3 mm continuum is basically unresolved. The mean flux 
density of 0.9 ±0.2 mjy measured from the integrated spec- 
tra (Fig. [H| is consistent with the 1.02 ± 0.07 mjy from a 
Gaussian fit to continuum map. 

Symmetric Gaussians directly fitted to the uv data give 
FWHMs of 0.85" and 1.35" for the projected intrinsic size 
of the HCN and HCO + lines respectively (we note that, 
as before, because these sizes are derived from the uv data, 
they are free of beam convolution). The molecular lines trace 
almost the same region as testified by the fits in spatial co- 
ordinates (Fig [7] central and right panels): the HCN emis- 
sion has an extent of 1.92" x 1.02" at PA=80°, similar to 
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Figure 8. HCN (upper) and HCO + (lower) integrated spectra of 
NGC 4051. Compared to NGC 2273 (Fig|2J and NGC 3227 (FigHJ 
both the continuum level and mean line width of FWHM=82 ± 
10 km/s are a factor of 2 lower. Lines are as in Figf3] 



1.65" x 1.27" at PA=72° for the HCO+ emission. Both the 
HCN and HCO + emission are clearly resolved when com- 
pared to the beam size. The velocity dispersion, obtained 
by fitting the spectra in Fig. [8] is about half that of either 
NGC 2273 or NGC 3227. The FWHM is indeed 74±10 km/s 
and 90 ± 10 km/s for the HCN and HCO+ lines respectively. 

In NGC 4051 the line fluxes are 0.91 ±0.05 Jy km/s for 
HCN, and 2.07 ± 0.05 Jy km/s for HCO+, both integrated 
over the central 3". This means that the HCN/HCO + inten- 
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Figure 9. The HCN emission map (rms = 
0.07 Jy beam^km s" 1 ) of NGC 6951. Labels are as in 
Fig. [JJ Contour levels are at 3, 4, and 5 times the noise level. 
The line emission appears at best only marginally resolved. 
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Figure 10. Integrated spectrum of NGC 6951 showing the HCN 
line (a fit for which is overdrawn in blue), with its large velocity 
dispersion (FWHM=173 km/s). 



sity ratio in NGC 4051 differs significantly from NGC 2273 
where the ratio is close to 1. This aspect of the line emission 
will be assessed, together with ratios of other line transi- 
tions, in a future paper. 



3.4 NGC 6951 

The properties of this source have already been presented 
by K07. As stated in that work, reconstructing an image 
with uniformly weighted visibilities provides high enough 
resolution that the emission from the ring does not con- 
taminate the nuclear signal. Here we briefly summarize the 
main characteristics of the nuclear emission that are required 
for our modelling, noting that the values listed in Table [3] 
agree very well with the K07 analysis and references therein. 
The HCN central emission (Fig. O originates from a re- 
gion of 1.44" x 1.10" at PA = 73°, at best only marginally 
resolved with respect to the beam (Table 1), with an in- 
tegrated flux of 1.02 ± 0.02 Jy km/s within 3". A single 
Gaussian is used to fit the line profile (Fig. I10[) . and yields 
a FWHM of 173 ± 29 km/s. 
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Figure 11. Channel maps of HCN and HCO + , superimposed on 
a greyscale image of the line. From top to bottom: NGC 2273, 
NGC 4051, NGC 3227 (left) and NGC 6951 (right). The red 
and blue channels are summed over different ranges from source 
to source, based on the line profile: ±(20 — 90) km/s (HCN) 
±(20 - 100) km/s (HCO+) for NGC 2273, ±(50 - 200) km/s for 
NGC 3227, ±(20 - 70) km/s (HCN) ±(20 - 80) km/s (HCO+) 
for NGC 4051, and ±(10 - 200) km/s for NGC 6951 (where the 
different velocity resolution of the data implies a different channel 
binning). The contour levels are at 2 a, 3<r, and 5 a for NGC 2273, 
NGC 3227 and NGC 6951, while for NGC 4051 are at 3a, 5 a 
and 7cr. The white 'plus' signs denote the peak of the continuum. 
In all cases the peaks of the red and blue channel maps are well 
separated, showing that the velocity gradient is resolved. 



4 PHYSICAL PROPERTIES OF THE DENSE 
GAS 

In this Section we first compare the dense gas kinematics to 
the warm gas and stellar kinematics. We then describe the 
dynamical models we apply to the data, under the assump- 
tion that the gas is distributed in a rotating disk. 

The high spectral resolution of the data enables us to 
measure a velocity gradient even though the sources are 
only marginally resolved, by generating red and blue channel 
maps - a technique that is equivalent to the spectroastrome- 
try used at optical wavelengths. It is a well-known technique 
that has been applied in many different circumstances, and 
is able to trace velocity gradients on spatial scales much 
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smaller than the spatial beam or see ing. The gain can eas - 
ily be a factor 10 or more (e.g. see iGnerucci et alj|20ld ). 
since the measurement simply depends on the precision with 
which the centroid of the emission can be measured, which 
is typically possible to ~ 1/10 of the beam size. The channel 
maps are generated by summing the line flux over positive 
and negative velocity ranges with respect to the line cen- 
ter (zero velocity). Figure [11] shows the channel maps of 
both HCN and HCO + for all targets and Table |3] lists the 
relative centroids separation and position angle of the vec- 
tor joining them. We note that the centers of red and blue 
channel maps are derived each time from the same dataset, 
and are spectrally close to each other. Therefore any offsets 
between red/blue emission due to phase calibration are neg- 
ligible and the uncertainties on their separation and relative 
position angle are only due to centroid positioning. For each 
object and emission line, we derive the errors on centroid 
position running Monte Carlo realizations of the red and 
blue channel maps, as described in Section 3. We then esti- 
mate uncertainties on their separation and relative position 
angle using a standard error propagation. In all the cases, 
the red and blue centroids are well separated, in most cases 
by a significant fraction of the beam sizes given in Table 1. 
The velocity gradients are thus well resolved. However, we 
note that the red and blue maps are not always centered 
symmetrically about the continuum, the location of which 
is indicated by the white plus signs on the maps (excluding 
NGC 6951, for which we have no continuum measurement). 
This is related to the apparent off-centre nature of the line 
maps in earlier figures, which may indicate that the HCN1- 
and HCO + 1-0 emission is not exactly axisymmetric, as 
might be expected if there are spatial variations in gas den- 
sity, temperature or excitation. 

For NGC 3227 the separation of the centres of the red- 
and blue-shifted emission is 0.38", and t he velocity gradient 
is oriented similarly to tha t of the stars (|Davies et al.ll2006l . 
120071 ; iBarbosa et alj I2OO6I . hereafter B06) and warm gas as 
traced by the H 2 (1-0) S(l) line (H09). Similarly, the cen- 
troid separation of 0.55" — 0.60" in NGC 4051 is consistent 
with the stellar kinematics (B06), but there is a mismatch 
in its position angle with respect to the warm gas and stellar 
kinematics in H09. This discrepancy can be understood in 
terms of both the angle measurement error (±10°), and the 
field of view of the IFU (0.56" x 2.24", H09) compared to 
the interferometer beam size in Table 1. For NGC 2273 we 
can compare the dense gas kinematics only with the stellar 
component (B06), and these are consistent. For NGC 6951 
stellar and warm gas kinematics are not available, and so we 
compare t he HCN kinematics to those of the ionized gas re- 
ported bv lStorchi-Bergmann etahl (|2007f ). The peaks of the 
red and blue HCN channel maps are separated by 0.60" at a 
PA consistent with the line of nodes inferred for the ionized 
gas. We can thus conclude that in our nuclei, the velocity 
gradients of the dense gas are, as far as can be ascertained, 
consistent with those observed with other tracers. 

The nuclear gas kinematics of our targets exhibit or- 
dered rotation. For example, high resolution observations of 
NGC 3227 and NGC 4051 show that the PA of the warm H 2 
remains constant throughout the nuclear region (H09). This 
suggests that the nuclear gas has no significant radial mo- 
tions (e.g., in/out flows) and no warp (i.e. there is no twist- 
ing of the kinematics major axis PA and/or inclination an- 




10 20 30 40 50 



Figure 12. HCN velocity dispersion map of NGC 6951, which 
shows a ridge of high dispersion along the minor axis (the typi- 
cal signature associated with severe beam smearing of a rotating 
disk), rather than along the radio jet (between the major and mi - 
nor axes) as has been reported bv lStorchi-B ergmann et al.l (120071 ) 
for the outflowing ionised gas. 



gle). The same conclusion is drawn in IBarbosa et ail (|2009h . 
where the sample includes also NGC 2273: although the ion- 
ized gas shows signs of outflowing motion, the kinematics are 
rather dominated by circular rot ation, coplanar with the 
stellar component. In NGC 6951, IStorchi-Bergmann et al.l 
l|2007h found evidence of an outflow due to the interaction 
of a radio jet f Saikia et al.l 120021 ) with the circumnuclear 
ionized gas. They measured a significantly enhanced veloc- 
ity dispersion (increasing from ~ 80 km/s to as much as 
150km/s) in two blobs of gas at PA 156°, an orientation 
close to that of the HCN velocity gradient. To verify if this 
feature is present in our HCN data, we show in Fig. [T5] the 
velocity dispersion map of NGC 6951. The observed disper- 
sion pattern has a ridge of high values perpendicular to the 
velocity gradient, which is the typical signature of a rotat- 
ing disk that is subject to severe beam smearing. As such, in 
NGC 6951, the radio jet appears to have negligible influence 
on the dense gas kinematics on these scales. We therefore 
conclude that there is no evidence for dense gas in radial 
motion in the sources considered here. Instead we expect 
HCN and HCO + to trace disk rotation. 

One intriguing issue is the differences between a mea- 
surements from different tracers. The velocity dispersion of 
the dense gas is systematically smaller than that of both 
the warm gas and stars. Once beam smearing is taken into 
account (see Sec 14. ip , this effect is even more pronounced 
and implies that HCN and HCO + have intrinsic velocity 
dispersions about 1/2 that of the warm gas and/or stellar 
component. This is perhaps to be expected, because of the 
high densities and lower temperatures traced by the HCN 
(HCO + ) molecules compared to the H2 1-0 S(l) line. One 
might speculate that there could be a stratified structure 
with the denser gas at lower scale heights (i.e. lower velocity 
dispersion). Alternatively, since the H2 (1-0) S(l) emission 
traces the edges of individual molecular clouds, if these are 
being ablated, then the line might preferentially trace out- 
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Input and Output Parameters of Dynamical Models 



Source 


Fixed 




Input 






Output 






e 


PA 


width 


height 




Major x Minor 


axis r/b sep 


FWHM 




deg 


deg 


// 


tt 


Arb. Uni. 


" X " 


// 


km/s 


NGC 2273 


51 


24 


1.5 


0.43 


2.9 


1.91 x 1.28 


0.99 


180 


NGC 3227 


55 


-37 


0.54 


0.14 


1.1 


1.27 x 1.05 


0.38 


207 


NGC 4051 


50 


29 


0.88 


0.25 


1.9 


1.88 x 1.09 


0.58 


88 


NGC 6951 


40 


-50 


0.99 


0.24 


2.1 


1.52 x 0.96 


0.62 


172 



Table 4. Dynamical simulations for a thick rotating disk. (1) Source name. The inclination of the molecular disk from H2 1-0 S(l) or 
stellar kinematics (2), and the PA (east of north) of the velocity gradient (3), are fixed parameters. The input parameters left free to 
vary are the disk diameter (4), thickness (5), and the scaling factor of the rotation curve (6). We adopt a Gaussian form for the radial 
mass distribution. The diameter and thickness are specified as FWHMs in, and perpendicular to, the disk plane. The model outputs are: 
the projected major and minor axes (7), the red and blue channel separation (8), and the line width (9). These can be compared directly 
with the observed quantities in columns (4), (6), and (8) of Table [3] For sources with two emission lines detected, the output values fall 
in between the observed quantities measured in Tab [3] for HCN and HCO + . 



flows. Since the ablated gas would form filamentary struc- 
tures, this would also tend to increase the velocity disper- 
sion. Although there is no sign of this in H09, it still remains 
a possibility. Despite the dispersion of the HCN (HCO + ) be- 
ing lower than that of other tracers, it still remains high. In 
order to understand this, we first assess how much of the 
observed line width is due to beam smearing of a velocity 
gradient. 

In the following we describe how we model the data, 
and how the model is constrained by four observed param- 
eters. When data for two molecular lines are available, we 
can simply take the mean values of the quantities listed in 
Table [3] Indeed, allowing for the various uncertainties (i.e. 
Monte Carlo realizations plus absolute position errors), all 
the quantities for the two lines in NGC 2273 and NGC 4051 
are almost similar. Once we have modelled the kinematics, 
we discuss the implications in the context of the AGN-star 
formation connection. 



4.1 Modelling the kinematics 

The aim here is to understand the basic intrinsic prop- 
erties of the emitting gas, specifically the origin of the 
linewidth. The model we use is the IDL code DYSMAL 
which is described in Appendix A of iDavies et al.l (120111 ) 
and has b een applied in a variety of cases in volving rotat- 
ing disks (IDavies. Tacconi. fc Genzel 2004al lbl; IDavies et al.l 
120091 : iGenzel et al.ll2008l ; ICresci et alTl2009h . 



DYSMAL is specifically designed to quantify the impact 
of the spectral and spatial beam smearing on an axisymmet- 
ric rotating disk, and thus allows us to infer the intrinsic 
properties of the dense gas from the observed ones. From a 
given set of input structural parameters, the code computes 
a data cube with two spatial and one velocity axis, which can 
be analysed in the same way as the original data. The ini- 
tial setup includes one or multiple mass/light components, 
parametrized by azimuthally symmetric functions (e.g., Ser- 
sic or Gaussian profiles; rings, etc.), the total rotationally 
supported mass, the scale height, inclination, and PA of the 
major axis on the sky. The model is then convolved with 
the beam shape and spectral resolution profile and sampled 
at the observed pixel scales. The total mass and its radial 



distribution fully define the rotation curve, which is com- 
puted assuming the disk is entirely supported by ordered 
rotational motion (i.e. a thin disk). An important param- 
eter for the dynamical model is the elliptical beam shape. 
This is crucial for modelling marginally resolved sources, in 
order to properly account for beam smearing. Interferomet- 
ric data, such those presented here, have a well characterized 
beam, which is listed in Table 1. 

DYSMAL can accommodate a departure from the ideal 
thin disk assumption, by allowing the thickness of the disk 
to be defined. If a finite scale height is specified, the associ- 
ated thick disk kinematics can be derived in 2 ways. With 
the first, the code makes an estimate of the local a based on 
the vertical height H and rotation velocity v at radius R, 
by setting R/H — v/cr as for a compact disk. Alternatively, 
once the geometry has been set, an additional fi xed and uni- 
form isotropic dispersion term can be included (jCresci et al.l 
120091 ; IDavies et aj]|201ll ) . 

We have four observable constraints related to the kine- 
matics: (i) the integrated linewidth; (ii) the major and (Hi) 
minor axis extent; and (iv) the separation of the red and blue 
channel map peaks. We can use these to match a dynami- 
cal model to the data via a x 2 minimization. The observed 
quantities in Table [3] can thus be compared with the output 
of our dynamical modelling, which is given in Tab [4] 

Due to the limited number of constraints, we need to 
make some reasonable simplifying assumptions as follows: 
(a) the gas exists in a rotating disk-like structure for which 
the inclination is the same as for the stars and/or for the 
H2 (1-0) S(l), and the position angle of which is set by the 
orientation of the line joining the centroids of the red and 
blue channel maps (PA r /f, in Tab [3] mean values are used if 
two lines are available) . These are given as the fixed param- 
eters in Table [4] (b) The mass and light follow the same dis- 
tribution. We use a Gaussian mass distribution, and derive 
the rotation curve from this. We note that assuming a uni- 
form or exponential mass distribution makes no significant 
difference to the results, (c) We can apply an axisymmetric 
model to emission that is not fully axisymmetric. In fact, 
significant velocity gradients due to unresolved non-circular 
motions can alter the kinema tic symmetries. As shown b y 
recent numerical simulations (|Hopkins fc Quataert lipoid ). 
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Toy models for NGC 3227 
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model 1 


0.5 


0.01 


0.5 


1.24 x 1.03 


0.50 


111 


model 2 


0.5 


0.01 


1.0 


1.24 x 1.03 


0.44 


191 


model 3 


0.5 


0.10 


0.5 


1.26 x 1.05 


0.37 


181 



Table 5. Toy models for NGC 3227 using thin and thick rotating disks. The inclination and position angle of the disk are 55° and —37° 
respectively. The input parameters, fixed in the three examples, are the disk diameter (2), thickness (3), and the scaling factor of the 
rotation curve (4). The shape of the rotation curve is set by adopting a Gaussian distribution for the mass profile. The diameter and 
thickness are specified as FWHMs in, and perpendicular to, the disk plane. The model outputs are: the projected major and minor axis 
(5), the separation of the centroids of the red and blue channel maps (6), and the molecular line width (7). 



the BH feeding process is related to gravitational instabili- 
ties that generically form lopsided, eccentric disks that prop- 
agate inwards from larger radii. This lopsided structures 
can dominate the di sk kinematics within the central 10 pc 
(Hop kins et al. ||2012| ). We consider it unrealistic to try and 
allow the model to account for details in the observed flux 
distribution. Instead, to overcome this limitation, we ignore 
the fact that the centroids of the red and blue channel maps 
are not symmetrically placed either side of the continuum 
peak, and simply use their relative position. 

To match a model to the observed data, several of the 
input parameters in Table [4] are left free to vary. These are 
the intrinsic width (ii) and height (H) of the rotating disk 
together with the scaling of the rotation curve. The last 
(M 3C aie) is defined by specifying, for some radius, an en- 
closed mass, defined as the mass that is supported solely by 
ordered rotation (and which may therefore differ from the 
actual mass for a thick disk). The accuracy of the centroid 
positions (and thus r/b separation) is, like the other uncer- 
tainties, estimated using Monte Carlo techniques. We create 
10 3 realizations of the model cube representing NGC 3227 
in Table 5, each time randomly adjusting the data values 
according to the measured S/N in Table [3l and re- fitting 
Gaussians. This yields ±0.02". We choose NGC 3227 for 
this test because it has the lowest S/N, and so for the other 
sources the uncertainties should be even smaller. 

To illustrate that the thin disk assumption is actually 
not appropriate to reproduce the observations, we show, in 
Table 5, several toy models for NGC 3227. The aim of these 
examples is not to fit the data (i.e. all the input parame- 
ters are not free to vary), but rather to demonstrate how 
the geometry and dynamics of the disk affect its observed 
properties. After smoothing with the spatial beam, the size 
in each case is comparable to that of the data, but there are 
large discrepancies with the channel map separation and the 
linewidth. We start in model 1 with a FWHM of 0.5". From 
a comparison with the NGC 3227 data in Table the size 
is about right, but the channel map separation is rather too 
large and the linewidth much too small. We therefore ad- 
dress in model 2 whether the discrepancies might be resolved 
by increasing the rotation rate. To do this, we increase the 
rotationally supported mass of the first model by a factor 
of 2. While this does broaden the line as required, it still 
does not match the r/b channel separation, which remains 
too large. Instead, in model 3, we put the mass scaling back 
to its initial value, and instead consider a thick disk. To do 
so requires an additional intrinsic dispersion which broadens 



the line width; and in addition the beam smearing effects on 
the thick disk reduce the r/b separation. The model yields 
a reasonable match for all the output parameters. 

Since a thick disk appears better suited to the data, 
we have left all 3 input parameters free to vary for each 
source. A comparison of the measurements in Table [3] with 
the model outputs in Table [4] clearly shows that the models 
can match all the observed constraints; and that the large 
velocity dispersion of the molecular gas is associated with 
its geometric distribution. Indeed, for all of our sources, the 
dense gas exists at significant scale heights above the disk 
plane. Table [S] lists the intrinsic quantities derived from 
models fitted to the dense gas kinematics. It shows that 
the ratio H/R (and thus <j/v) is about 0.3, if we allow DYS- 
MAL to estimate the dispersion associated with a given scale 
height. On the other hand, if we use a fixed uniform addi- 
tional a term, we obtain a/v ~ 0.5. These two ratios are in 
approximate quantitative agreement, given the contrast be- 
tween the simplicity of the models and greater complexity of 
the data. We thus adopt a mean a/v = 0.4 in the following 
discussion. 

The dynamical models allow us to compare the intrin- 
sic dense gas kinematics with those observed for the warm 
molecular hydrogen and stars. The intrinsic velocity disper- 
sions in Table [6] for NGC 3227 and NGC 4051 are about 
half of those measured from the warm gas, as given in Tab 3 
of H09. The same conclusion is reached by comparison to 
the stellar dispersion in Table 4 of B06. As mentioned in 
Section 4 the difference can be reasonably ascribed to the 
high densities traced by HCN and HCO + . Nevertheless, the 
dense gas is still scattered to significant scale heights. 

We note that it is also rotating faster than either the 
warm H2 or stars. Indeed the rotation velocity in our models 
at a radius of 30 pc is 105 km/s in NGC 3227 and 47 km/s 
NGC 4051, compared with 59 km/s and 37 km/s, respec- 
tively, measured from the H2 (1-0) S(l) line (H09). 

4.2 Comparison with other galaxies 

In this section, we briefly compare the intrinsic dispersion 
of the dense gas in these four AGN as given in Table [(J with 
that observed in two other extreme regimes. 

The center of the Milky Way (MW) is a good choice 
for a reference of a quiescent galaxy because of its proxim- 
ity. Thus, even with modest angular resolution and sensi- 
tivity, one can detect HCN on small physical scales and at 
low gas masses; and HCN emission has been targeted nu- 
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22 


6 
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NGC 6951 
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11 
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Table 6. Intrinsic kinematics of the modelled disk. (1) Source 
name. (2)-(3) Scaling radius and height given in parsecs. (4)-(5) 
Rotation velocity and velocity dispersion at R. 



merous times. IChristopher et all ((2005) observed HCN(l-O) 
and HCO + (1-0) in the circumnuclear disk of the MW and re- 
solved more than 20 molecular cores. The mean velocity dis- 
persion they measure is < a >= 11 ±0.4 km/s, smaller than 
we find in our AGN. However, the MW circumnuclear disk 
is only 3 pc across, rather less than the size scales our data 
are p robing. On more comparable scales of ~ 100 pc in the 
MW, |Lei (|l996t ) found typical values of a ~ 10 - 21 km/s, 
which is still smaller than in three of the four AGN. This 
suggests that there may be a heating source in the central 
regions of AGN that maintains the high dispersion in those 
objects, and which is not present in the MW. 

The opposite regime is represented by strong star form- 
ing galaxies also harbouring AGN, such as Ultraluminous 
In frared Galaxies (ULIRGs). In the thre e ULIRGs observed 
bv llmanishi. Nakanishi. fc Kohnd (|2006t ). the HCN(l-O) and 
HCO + (1-0) lines show extremely large line widths, in the 
range 140-300 km/s FWHM. These are the directly observed 
integrated values, which should therefore be compared with 
our 80-200 km/s from Table [3] The similar ranges sug- 
gest that the central regions of AGN may be as extreme 
as ULIRGs. However, we note that the beam smearing cor- 
rection is unknown for the ULIRG measurements, and since 
the data trace large spatial scales (~ 1 — 3 kpc), this may 
partially explain the high velocity dispersions observed. 

The implication of discrepant velocity dispersions in dif- 
ferent regimes is intriguing, as it might reveal a link between 
the central disk structure and the nuclear star formation. 
This issue is discussed in the next section. 



5 DISCUSSION 

We have found that in our four active galaxies the intrin- 
sic HCN and HCO + velocity dispersion lies in the range 
~ 20 — 40 km/s. While this is about half of the velocity 
dispersion measured for the warm H2 (H09), it is still re- 
markably high. Indeed, the high velocity dispersion of the 
gas implies that the dense gas distribution is geometrically 
thick. But the mechanism responsible for maintaining it is 
unclear. One inevitable conclusion is that such structures 
have to be clumpy, and the high dispersion of the clouds 
must be supported by a force other than thermal pressure, 
otherwise inelastic collisions should collapse the structure 
to a thin disk within an orbital timescale (a few Myr at 
our sc ale radii). Supernovae explosions in a p ost-starburst 
phase (|Wada. Papadopoulos. fc Spaansl 120091 ) appear able 
to explain the lower end of the velocity dispersion range 



Table 7. Gas properties within the scaling radius R. (1) Dy- 
namical mass M^yn = (v 2 + 3cr 2 )R/G (see text for details). Ex- 
cluding the uncertainty of the coefficient for the a term, ran- 
dom errors on dynamical masses are 15 — 20%. (2) HCN lumi- 
nosity, which h as been found to be proportional to t he dense 
gas mass (EqrTl lGao fc Solomodl2004lKrips et al.ll200Sl 1. (3) Gas 
fraction, for which the relative error is in the range 20 — 30%. 

4) Toomre Q parameter norm alized to the critical value of 0.68 

Dekel. Sari, fc Ceverin"o1l2009l') . 



(<r ~ 20 km/s). On the other hand, the higher disper- 
sions require an alternative process. Gravitational turbu- 
lence in duced by external gas accretion into the nuclear 
region (|Vollmer. Beckert. fc Daviesj |200S| . hereafter V08), 
could represent a plausible way to explain values of 40 km/s, 
if this mechanism can be shown to work. Three phases char- 
acterize the V08 evolutionary model: first a short and mas- 
sive gas inflow from the host generates a massive turbulent 
disk where star formation occurs. In a second phase super- 
novae explosion clear the intra-cloud medium leaving a ge- 
ometrically thick disk dominated by dense gas cloud cores. 
Finally, when the mass accretion rate into the central regions 
decreases sufficiently, the disk becomes thin and transparent. 
To assess the evolutionary phase characterizing our sources, 
we first estimate whether or not the modelled disks are sta- 
ble against star formation, and we then investigate the SF 
properties in their nuclear regions. 



5.1 Can stars form in the dense gas? 

A rotating gaseous disk becomes unstable once the local 
gravity overcomes both differential rotation and turbulen t 
pressure, i.e. when the Toomre Q parameter (|Toomrelll964l ) 
is smaller than a critical value Q c : 



tvGT, 



< Qc 



(1) 



In Equation U G is the gravitational constant, a the lo- 
cal velocity dispersion; and k is the epicyclic frequency, 
which is a function of the angular velocity fl. An isother- 
mal thick disk has a critical Q parameter Q c = 0.68 
(IDekel. Sari, fc Ceverinol I2009T ). For a uniform disk the 
epicyclic frequency becomes: 



= VSv/R, 



(2) 



with v the rotation velocity at radius R. The surface density 
of the gaseous disk in Eq. [T]is: 



Sgas — 



Mqas 



fg a s -Mdy n 



(3) 



7T.R 2 7T.R 2 ' 

where, in the right part, we replace the gas mass M gas with 
the derived dynamical mass Md yn scaled by the gas frac- 
tion fgas - Because of significant random motions (a/v ~ 0.4 
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from Table [6]), we note that Mdyn, in turn, depends on both 
rotation velocity and dispersion: 



M dyn = {v 2 + 3a 2 )R/G. 



(4) 



Equation [4] arises from combining the kinetic energy of the 
random motions, which we take to be 3<r 2 (assuming a comes 
from macroscopic turbulence, and that we measure a one- 
dimensional component of an isotropic three-dimensional 
distribution), with a v 2 term for the rotation. Finally, we 
use the relation between size scales and kinematics for a 
compact disk, which gives: 



vH/R. 



(5) 



By substituting Equations [21 EI SI an d into the Q param- 
eter definition (Eq. [JJ we obtain a simple expression for Q, 
that depends only on the geometric structure of the thick 
disk and the gas fraction: 



Q = \/3 



H 



'R(l + 3H2/R*)f gaa - (6) 

We note that Eq.[6] is, for example , analogous to the expres- 
sion derived in lGenzel et al.l ([2011), but with an added term 
corresponding to the mass supported by random motions. 

It is straightforward to obtain the dense gas fraction 
fgas = M g as/Mdyn, once the molecular gas mass is esti- 
mated. The HCN luminosity (Lhcn) is related to M gas 
through l|Gao fc Solomon! |2004 iKrips et all l200Sh : 



Mg aa = WLhcn Mq (Kkms pc 



(7) 

The dynamical mass, as given in Eq. [4] is listed in Ta- 
ble [7] together with the HCN luminosity, gas fraction and 
the Toomre Q parameter computed at the scale radius R 
taken from Table [6] The estimate of the uncertainties given 
in the Table is based only on the propagation of random 
errors, and does not include systematic uncertainties such 
as the coefficients used to estimate the dynamical mass or 
the uniqueness of the relation between HCN luminosity and 
gas mass. Nevertheless, with a gas fraction of ~ 10 — 20% as 
derived above, there are indications that the Q parameter 
is greater than the critical value, which would imply that 
the central 30-100 pc of our AGN are stable against star 
formation. 



5.2 Is there evidence for star formation? 

The high Toomre Q parameters derived for our sample 
(Table [7]), indicate that it should not be possible for the 
gas to rapidly form stars. We expect, therefore, that no 
significant star formation should be detected in the cen- 
tral ~ 100 pc regions observed with the PdBI. To ver- 
ify this conjecture, we look for ongoing SF by investi- 
gating mid-infrared spectra in the context of other obser- 
vations. The 5—15 fim waveband is expected to show 
strong polycyclic aromatic hydrocarbon (PAH) features if 
the m olecules are excited by radiation from young massive 
stars (|Allamandola. Tielens. fc Barked 1 19891). The mi d-IR 



spectra obtained with Spitzer /IRS i Houck et al.|[2004T ) are 
shown in Fig. 1131 All the sources exhibit strong PAH emis- 
sions, suggestive of significant star formation. But it is not 
clear where the SF is spatially located. Indeed, due to the 
IRS slit aperture of 3.6", the spectra are integrated over the 
central 450 pc, 305 pc, 180 pc and 370 pc respectively for 



NGC 2273, NGC 3227, NGC 4051 and NGC 6951. Data ob- 
tained at high angular resolution can help, but it can still 
sometimes be hard to reach a firm conclusion. 

As described in Section 3.2, in NGC 3227 the HCN 
nuclear emission is surrounded by an annular structure 
seen wit h CO (2-1) line and H-band continuum observations 
(Fig. [5] ISchinnerer. Eckart. fc Tacconil l2000al . iDavies et al.l 
120061 ) . This circumnuclear ring hosts the bulk of the current 
SF; while at radii < 30 pc, SF has recently cea sed. This has 
been shown by adaptive optics observations (jDavies et al.1 
2006), which reveal that the Br7 equivalent width is incon- 
sistent with on-going SF. Based on this and other diagnos- 
tics, these authors concluded that the inner tens of parsecs 
in NGC 3227 have experienced a short starburst event in 
the last 40 Myr, which has now ceased. 

For both this source and NGC 2273, B06 find a re- 
duction in the stellar velocity dispersion across the ring, 
which they ascribe to recent SF in dynamically cool gas. 
The implication is that the higher dispersion in the nucleus 
indicates that SF is not active there. For NGC 2273, the 
circumn uclear ring is clearl y visible also in HST/NICMOS 
images dMartini et all 20031). From t heir spectroscopic data, 
iFerruit. Wilson, fc Mulchaevl (|2000l ) relate Ha emission to 
HII regions in the ring, whereas they argue that the nu- 
clear [OIII]/Hq ratio is typical of AGN. Thus, although 
we cannot firmly rule it out , these results from B06 and 
IFerruit. Wilson, fc Mulchaevl (|2000l ) suggest that, as for 
NGC 3227, there is no current SF in the nucleus. 

In NGC 4051 stellar kinematic measuremen ts have been 
preve nted by the strong AGN emission (B06; Riffel et al.l 
120081 . hereafter R08), and the central physics are more puz- 
zling. The Br7 emission is very compact and it is hard to 
draw firm conclusion: R08 find no rotation and argue that 
most of the emitting g as is not confined to t he pl ane of 
the galaxy; in contrast, iMiiller-Sanchez et all (|201 lh trace 
a Br7 velocity field consistent with rotation, that is seen 
also in the H2 1-0 S(l) and stellar kinematics (H09). From 
the H2 line ratios R08 argue that the gas is thermalised 
rather than fluorescently excited as one might expect from 
recent star formation. Although this could instead be due 
to high densities, the H2 (1-0) S(l)/Br7 ratio in the nucleus 
is about 1, higher than expected if the heating process is 
due to UV radiation from young stars. Moreover, R08 also 
detect [OIII] and [CaVIII] lines characterized by blue wings 
and argue that significant line emission originates from the 
AGN narrow-line region. Thus, a similar origin of the com- 
pact Br7 emission is plausible. On the other hand, if we do 
attribute the Br7 emission to SF, we c an estimate its equiva- 
lent w idth once the flux measured bv lMuller-Sanchez et all 
|201ll ) is corrected for dilution (see Fig. 1 in R08). This 
yields an equivalen t width of ab o ut 7A , a low value for SF 
that, according to Davi es et al. J2007f). implies ceased SF 
activity. Finally, iRodrfguez-Ardila fc Viegasf (|2003u are able 
to give only an upper limit for the 3.3 /im PAH feature using 
a 0.8" slit. We can thus conclude that the SF is NGC 4051 
is mainly located between 0.8" and 3.6", and if a nuclear 
component within the central 40 pc is present, it is not sig- 
nificant. 

In addition to its AGN, NGC 6951 also exhibits a 
pronounced starburst ring located at 5" (480 pc) from 
the nucleus. Plausibly, most of the PAH emission may 
originate from this ring, which has been observed in 



© 2012 RAS, MNRAS 000.[THT5l 



14 E. Sani et al. 



1.00 



0.10 



0.01 



NGC2273 
NGC3227 

NGC4051 



NGC6951 




10 

Rest wavelength [/j.m] 



15 



Figure 13. Spitzer /IRS spectra for the four AGN. The spec tra are arbi t rarily scaled for clarity, and are ordered with decreasing AGN 
luminosity from the top. The data reduction is described in I Sani et al.l d201GT) ■ The fo ur spectra show typical PAH em ission features 
at 6.2 )im, 7.7 [ im, 8 .6 (im, 11.3 ^ m, 12.7 )im and have been previou sly analysed by |Diamond-Stan ic fc Riekel J2010f) (NGC 2273), 
iMelendez et al] l|200Sh (NGC 3227). ISani et al.l lj20ld l (NGC 4051), and lDiamond-Stanic fc Riekel (201(1 ) (NGC 69511. 



several wavebands: near- IR jMartini et all 120031 ). optical 
(e.g. Istorchi-Bergmann et alj 20071 and references therein), 
and radio ( Saikia et al.l I2002T ). High resolution observa- 
tions detect strong CO and H CN emission both in the 
ring and the central AGN (K07; iGarcfa-Burillo et al.ll2005l. 



iKohno. Kawabe 



Vila-Vilard 1 19991 ). While the HCN /CO 
ratio in the ring is typical of a starburst, the HCN is 
significantly enhanced in the nucleus. The higher nuclear 
ratio could be due either to denser and/or hotter gas 
than in the starburst ring, or because the gas chemistry 
in the nucleus of NGC 6951 is dominated by X-ray ra- 
diation from the AGN, skewing relative abundances (see 
iDavies. Mark, fc Sternberg! 12012). Thus, again, significant 
SF close around the AGN cannot easily be excluded. The 
radio observations can help to address this issue. The super- 
nova rate, esti mated from the n uclear radio flux density, is 
~ 0.003 yr _1 |Saikia et al.ll200"2T ). much smaller than would 
be expected for strong ongoing SF, and rather less tha n in 
NGC 3227 where the SF has ceased (|Davies et al.ll2006l ). In- 
deed, if the SF is continuous, we ca n estimat e an a ssociated 
SF rate (SFR) from Equation 20 in lCondonl (| 19921 ). obtain- 
ing SFR(M> 1M©)~ 0.15 M yr _1 (including sta rs below 
5 Mq , thus a factor of 2 higher as for the original ICondonl 
(1992) relation). This is comparable to the t ime-averaged 
SFR for the nuclear regions of other AGN (D avies et al.l 
120071 ). and provides an upper limit to the actual on-going 
SFR. 



6 CONCLUSIONS 

We have analysed 3 mm interferometric data obtained with 
the PdBI tracing the HCN(l-O) and HCO+(1-0) molecular 



lines in four local AGN. Our main conclusion, that these 
lines typically trace a thick disk with suppressed star for- 
mation, is based on the following: 

1. The line emission is marginally resolved in all four 
Seyfert galaxies. When data for both HCN and HCO + are 
available, their spatial extents are in good agreement within 
the errors. 

2. The global kinematics are consistent with rotation. 
We create dynamical models that match the observations, 
using four constraints, namely the integrated linewidth, the 
diameters of the major and minor axes, and the separation 
of the red and blue channel map centroids. The resulting 
structure is a geometrically thick disk. 

3. The most remarkable feature is a high velocity dis- 
persion that, once corrected for beam smearing using the 
dynamical models above, lies in the range 20-40 km/s. This 
is about a factor 2 lower than the stellar and/or H2 1-0 S(l) 
dispersion, but is still remarkably high considering that we 
are tracing dense molecular gas. It implies a vertical height 
H~ 30% of the scale radius R. 

4. Finally, we estimate the gas fraction within R (f ga s ~ 
10 — 30%) and the Toomre Q parameter, finding that the 
central ~ 100 pc should be stable against star formation. 
This result appears to be observationally confirmed by the 
lack of ongoing nuclear SF in the nucleus (although there is 
plentiful evidence for on-going star formation in circumnu- 
clear rings). NGC 4051 remains the most puzzling source, 
with the lowest velocity dispersion and unclear evidence for 
or against nuclear star formation. 
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